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1 Introduction 



The semiclassical methods in solving the Schrodinger problem are of extreme 
importance in understanding the global behaviour of eigenfunctions and en- 
ergy spectra, especially as a function of some external parameter, since usu- 
ally they are the only approximation known in the form of an explicit formula. 

The leading semiclassical approximation is just the first term of a certain 
^-expansion. The method goes back to the early days of quantum mechanics 
and was developed by Bohr and Sommerfeld for one-freedom systems and 
separable systems, it was then generalized for integrable (but not necessarily 
separable) systems by Einstein (1917), which is called EBK or torus quan- 
tization. In fact, Einstein's result was corrected for the phase changes due 
to caustics by Maslov (1961; see also Maslov and Fedoriuk 1981), but the 
torus quantization formula thus obtained is still just a first term in a certain 
^-expansion, whose higher terms are unknown in systems with more than one 
degree of freedom. Thus recently it has been observed (Prosen and Robnik 
1993, Graffi, Manfredi and Salasnich 1994) that this leading-order semiclassi- 
cal approximations generally fail to predict the individual energy levels (and 
the eigenstates) within a vanishing fraction of the mean energy level spacing. 
This conclusion is believed to be correct not only for the torus quantization 
of the integrable systems, but also in applying the Gutzwiller trace formula 
(Gutzwiller 1990) to general systems, including the completely chaotic ones, 
c.f. Gaspard and Alonso (1993). Therefore a systematic study of the ac- 
curacy of semiclassical approximations is very important, especially in the 
context of quantum chaos (Casati and Chirikov 1995, Gutzwiller 1990). This 
to end in full generality is an almost impossible task, but in some special 
cases it is possible to work out the quantum corrections to higher or even 
all orders (Degli Esposti, Graffi and Herczynski 1991, Graffi and Paul 1987, 
Salasnich and Robnik 1996, Robnik 1984, Narimanov 1995). On the other 
hand, in systems with one degree of freedom a systematic WKB expansion 
is possible at least in principle, and in a few cases can be worked out even 
explicitly to all orders, resulting in a convergent series whose sum is identical 
to the exact spectrum (Dunham 1932, Bender, Olaussen and Wang 1977, 
Voros 1993, Robnik and Salasnich 1996). 

Our goal in the present paper is to deal systematically with the WKB 
expansions for the angular momentum problem and for the Kepler problem. 
This is important not only from the point of view of mathematical physics 
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(formal existence of the systematic series, its convergence properties and 
the summation), but also because the Kepler problem is so fundamental in 
physics. To the best of our knowledge a detailed analysis of this problem has 
not been undertaken in the literature so far. 

We shall work out some next to the leading terms for the Kepler problem 
and show - under a conjecture about the higher terms - that exact result is 
obtained after all corrections have been taken into account and the result- 
ing series has been summed. This is nontrivial, because we know that the 
torus quantization of the 3-dim Kepler problem yields exact result, whereas 
the individual torus quanization of the radial and of the angular momentum 
problems is not exact. Thus our present work is the first systematic semi- 
classical expansion of the angular momentum problem as a pre-requisite to 
the full study of the 3-dim Kepler problem. 

To define the language and to introduce the notation we first give the 
essential formulas of the torus quantization. The Hamiltonian of the 3-dim 
Kepler problem is given by 

2 2r 2 r K J 

where 
and 

P* = L z (3) 

are constants of motion. Of course, the Hamiltonian is a constant of motion, 
whose value is equal to the total energy E. 

It is well known that the exact energy spectrum can be obtained with the 
Bohr-Sommerfeld (torus) quantization. To perform the torus quantization it 
is necessary to introduce the action variables 

h = ^jP^<P = P^ (4) 

I e = ^fPedO = L-I 4> , (5) 
I r = — I p r dr = - L. (6) 
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The Hamiltonian as a function of the actions reads 



-a 2 



H 2[J r + J fl + (7) 



and after the torus quantization 



1 1 

I r = (n r + -)k, I e = (n e + -)h, 1$ = n^k, (8) 



the energy spectrum is given by 



-a 2 



Enrl 2h 2 [n r + l + ir (9) 

where I = rig + n<p. (Each of the three quantum numbers is a nonnegative 
integer, and so is I.) This is the exact energy spectrum, which can also be 
obtained by solving the Schrodinger equation. Note that we have quantized 
the angular momentum L = I r + I e with a semiclassical formula L — (I + 
l/2)h. If we use the exact quantization of the angular momentum, i.e. L = 
+ 1), we obtain a wrong formula. How to explain this observation? 
In section 2 we treat the angular momentum problem by calculating the 
corrections to the leading torus quantization term, and in the section 3 we 
then proceed with the analysis of the radial Kepler problem, again by calcu- 
lating the corrections to the leading torus quantization term, now using the 
exact result for the quantized angular momentum. In section 4 we discuss 
the results and draw some general conclusions. 

2 WKB expansion for the angular momen- 
tum 

We consider the eigenvalue equation of the angular momentum 

L 2 Y(e,<f)) = \ 2 h 2 Y(e,<p), (io) 

where L 2 is formally given by the equation (2) with 

p ^ = - h ^W 2+coH6) ¥e h 
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p i- ti W (12) 

We can write the eigenfunction as 

Y(9,<f>) = T(9)e in ^, (13) 

and we obtain 

P*Y(0,<f>)=nlh 2 Y(9, ( j ) ), (14) 

and also 

T"{9) + cot (9)T'(9) + (A 2 - = 0. (15) 

Notice that h does not appear in this equation anymore. To perform the 
WKB expansion we introduce a small parameter e, which might be thought 
as proportional to H, and consider the eigenvalue problem 

e 2 T"(0) + e 2 cot (9)T'(9) = Q(9)T(9), (16) 

where 

n ^ 

Q{9) = W{9) - A 2 = - A 2 . (17) 

Thus small e limit is equivalent to the large and/or large A limit. The pa- 
rameter e helps to organize the WKB series; we set e = 1 when the calculation 
is completed. First we put 

T(0)=exp{is(0)}, (18) 

where S(9) is a complex function that satisfies the differential equation 

S' 2 {9) + eS"{9) + e cot (9)S'(9) = Q{9). (19) 
The WKB expansion for the function S{9) is given by 

oo 

S(9) = J2^ k Sk(9), (20) 

k=0 

and by comparing like powers of e we obtain a recursion formula (n > 0) 

S'l = Q, E S' k S' n . k + 5" n _! + cot (0)S" n _i = 0. (21) 

fc=0 



5 



Straightforward calculations give for the first few terms 

S'o = ~QK (22) 
^ = _ I cot (9), (23) 

S> 2 = -lQ^Q-5/2_ 1^(^-3/2) .1^2(^-1/2 

- \(± cot (9))Q~^. (24) 
The exact quantization of the wave function (18) is given by 

/oo 
dS = ^2 j> dSk = 2ni rig, (25) 

k=0 

where we have now set e = 1. This integral is a complex contour integral 
which encircles the two turning points on the real axis. Obviously, it is 
derived from the requirement of the uniqueness of the complex wave function 
T (Dunham 1932, Bender, Olaussen and Wang 1977). 
The zero order term is given by 

j> dS = 2ij dr^X 2 - W(6) = 2m(X - n ), (26) 

and the first term reads 

<j> dSi = ~ lnQ\ contour = -iri. (27) 

Evaluating InQ once around the contour gives 4ni because the contour en- 
circles two simple zeros of Q. 

All the other odd terms vanish when integrated along the closed contour 
because they are exact differentials (Bender, Olaussen and Wang 1977). So 
the quantization condition (25) can be written as 

oo 1 

J2fdS 2k = 2Tri(ng + -), (28) 

fc=0 
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and thus it is a sum over even-numbered terms only. The next non-zero 
term is given by 

= w ' 2[e) + \° I 



ud(\YJ ft* - w{0) - \v(0) 

+ l/*_j£(2_]. (29, 

These three integrals give (see the Appendix A) 

fdS 2 = ^, (30) 

where, importantly, the dependence drops out now. Thus up to the second 
order in e the quantization condition reads 

A + ^=/ + i (31) 

where / = rig + n^. The term 1/8A is the first quantum correction to the 
the quantization of the angular momentum. From this result we can argue 
("conjecture by educated guess") that the e 2k term in the WKB series is 
(k > 0) 

' dS 2k = 2m(f\2- 2k \ 1 ~ 2 \ (32) 



so that the WKB expansion of the angular momentum to all orders is given 
by 

£(!) 2_2fcAi " 2fc = z+ i ^ 

This is the exact formula for the relationship between / and A, because 

g (\\ 2 -2 k) l-2k = I V / 1+4A 2 ; (34) 

fc=o W 2 

and the equation y/l + 4A 2 /2 = 1 + 1/2 can be inverted and gives A = 
+ This completes our investigation of the semiclassical expansion for 
the angular momentum, where it remains in general to prove the conjectured 
formula (32) for k > 2. 
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3 WKB expansion for the radial Kepler prob- 
lem 

We consider the Schrodinger equation for the radial problem 

' h " °" +V{r)\i){r) = Ei){r) (35) 



L 2 dr 2 

where 



V{r) = ^ " f (36) 
We can always write the wave function as 

i 

ip{r) = exp{-<x(r)}, (37) 

where the phase a(r) is a complex function that satisfies the differential 
equation 

*' 2 (r) + £y>(r) = 2(E-V(r)). (38) 

t 

The WKB expansion for the phase is 

oo *■ 

o(r) = £(-)V fc (r). (39) 
k=o 1 

Substituting (39) into (38) and comparing like powers of h gives the recursion 
relation (n > 0) 

n 

a'l = 2(E - V(r)), £ a' k a' n - k + a" n ^ = 0. (40) 

k=0 

The quantization condition is obtained by requiring the uniqueness of the 
wave function 

/°° h r 
da = ]T(-) fc f da k = 2nn r h (41) 
fc=0 1 J 

where n r > 0, an integer number, is the radial quantum number. 

The zero order term, which gives the Bohr-Sommerfeld formula (6), is 
given by 

jdav = 2 J dr^/2(E -V(r)), (42) 
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and the first odd term in the series gives the Maslov corrections (Maslov 
index is equal to 2) 

j> da ± = -irh. (43) 

All the other odd terms vanish when integrated along the closed contour 
because they are exact differentials (Bender, Olaussen and Wang 1977). So 
the quantization condition (41) can be written 



00 h r 1 

E(-) 2fc f da 2k = 2n(nr + -)h, (44) 
k=o 1 ' z 



thus again a sum over even-numbered terms only. The next two non-zero 
terms are (Bender, Olaussen and Wang 1977) 

■ h vlj *2 1 d 2 f, V' 2 (r) 

- * da 2 = ~n 2 — -— / dr =, 45 

12 dE 2 J \J2(E - V(r)) 

^4i,. _ fc 4 r l & V"\r) 1 & [ dr V'\r)V"{r) 1 



(") 4 f da, = h[—— 3 J dr-j=±—- — m j J 2{E _ v{r)) i 

(46) 

A straightforward calculation of these terms gives (see the Appendix B) 
and 

Up to the fourth order in h the quantization condition reads 

So we have obtained the first two quantum corrections to the torus quanti- 
zation of the radial Kepler problem. Obviously at this point of truncating 
the series we get wrong spectrum if we use the torus quantized angular mo- 
mentum L = (/ + l/2)h, and this is still true if the series is expanded to all 
orders. However, for the anticipated infinite series expansion we shall obtain 
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the exact quantized value of the eigenenergies when using the exact angular 
momentum L 2 = 1(1 + l)Ti 2 . To show this we note that higher-order correc- 
tions quickly increase in complexity but each integral gives a polynomial in 
E with leading term E M , where M is the power of the operator d M /dE M 
in front of the integral (Barclay 1993). Differentiating M times leaves a con- 
stant independent of E. Since this happens in all terms in the series (with 
k > 0), the WKB corrections to the Bohr-Sommerfeld formula have no In- 
dependence. From this result we can guess the general formula, based on our 
two correcting terms to the torus quantization, namely 

rv 1 00 /!\ 1 00 /I\ 



where A = L/h, and so the % 2k term in the WKB series is (k > 0) 



(50) 



{^) 2k jda 2k = -2^^2- 2fe A 1 - 2fc . (51) 

In conclusion, the energy spectrum of the WKB algorithm to all orders is 
given by 

e wkb = Z^L (52) 

HrX 2^[(n r + |)+E^o(|)2- 2fe A^] 2 - 

Now, by using the formula (33) of the WKB expansion of the angular momen- 
tum, we obtain indeed the exact result E^^ B = E nr i, as given in equation 
(9). 

We can summarize the mathematical reason for exactness of the torus 
quantization formula (derived in the Introduction) for the 3-dim Kepler prob- 
lem: Since the problem is separable, the wave functions (for the angular 
momentum and for the radial part) multiply and their phases have the ad- 
ditivity property, and therefore the total phase written as | (a — ihS) must 
obey the quantization condition (uniqueness of the wave function). From 
the two formulae (32) and (51) one can see that the quantum corrections 
(i.e. terms higher than the torus quantization terms) do indeed compensate 
mutually term by term. 
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4 Discussion and conclusions 



In the present paper we offer (to the best of our knowledge) the first calcula- 
tion of the higher WKB terms beyond the torus quantization leading terms 
for the angular momentum and the radial Kepler problem. This analysis 
explains the curious compensation of the higher order quantum corrections 
(of the two separated problems) resulting in the exactness of the torus quan- 
tization for the entire 3-dim Kepler problem (see the Introduction). We have 
no reason to doubt that our conjectured general formulae (32) and (51) are 
correct for all k > 0, but this still has to be proved. 

We consider this kind of studies as important in understanding the accu- 
racy of the semiclassical methods, and much of the results in this context for 
1-dim problems are known, including some more general families of 1-dim 
potentials studied by Barclay (1993) which are characterized by the factor- 
ization property (Infeld and Hull 1957, Green 1965). One important future 
project is to analyze a more general class of the 1-dim potentials and in par- 
ticular to extend results to systems with two or more degrees of freedom, 
even if they are integrable (but not separable). Further, it remains as an 
important project to assess the accuracy of much more general (although 
mathematically not yet completely satisfactory, due to the divergent series 
expansions) methods like the Gutzwiller theory (1967-71, 1990), applicable 
to nonintegrable systems, including the chaotic systems (Gaspard and Alonso 
1993). 
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Appendix A 



In this appendix we show how to obtain the formula (30). In all integrals 
of this section the limits of integration are between the two turning points. 
After substitution z = tan (9), we have 

W"(6) _ K f dz (l + z>) _ 



fa . w m = m * f. 

J . /\2 _ W(f)\ . \ 2 - w 2 J 



- W{d) ^A 2 - ny ^ \z 2 -p 

' T (A 2 -n ) 2 + 2vrn (A 2 -nJ), (53) 



2n, 

where (5 = n|/(A 2 — nV), so that 



92 r M = * (54) 



d(\ 2 ) 2 J ^/A 2 - W(6) V 
For the other integrals we use the same procedure. 



from which we obtain 

9 f d0 ^m * (56) 



«9(A 2 ) y ^A 2 - W{9) n <t> 
The last integral gives 



r Jn cot2 ( ) 1 /" 7 1 / Z 2 , 1 1. . ^ 

/ d6 — = = . / dz—- — -z = tt( -). (57) 

J - ^(0) Jx^4 J * 2 (i + * 2 )V* 2 -/? x y 

In conclusion 

r 1 37T 1 . 7T . 1/1 1 N , 7ri . 

^2 = -?—— + - + T =7T- 58 

L 12n 2 V n/ 4 A 4A 
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Appendix B 



In this appendix we show how to obtain the formulas (47) and (48). In this 
section again all integrals are taken between the two turning points. For the 
first one, after substitution y — 1/r, we have 

[ ir v ' 2{r) = I ^v-^w+qy m 

J p{E - V(r)) J Wa + by-y* ' 1 ' 

where a = 2E/L 2 and b = lajl? . We observe that 

4 

h = I dy V = T^(48a 2 + 128ab 2 + 356 4 ). (62) 

Because we must apply the operator d 2 /dE 2 and a = 2E/L 2 , the only non- 
zero contribution stems from the integral 1^ and we obtain 



dE2 J \f%E - V(r)) L 
In conclusion we have 

To obtain the formula (48) we proceed in the same way. 

f dr V"\r) _ f d 9LV - 12LW + 4a V 
i Jo(R- V(r\\ J LJa + by - y' 2 



its leading integral is 

J 6 = / dy — V% = = -^(320a 3 + 1680a 2 6 2 + 1260a6 2 + 2316 6 ), (66) 
J v a + by — y 2 1024 
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from which we obtain 



& 3 V"\r) _135tt 



dE 3 J ^2(E-V{r)) ' L 3 ' 
For the last integral we have 

, V'\r)V"(r) r , 3LV - 8L 4 y 7 + 7L 2 a 2 y 6 - 2a 3 y 5 
dr i = / dy- 



y/2(E-V(r)) J Wa + by - y 2 

its leading integral is 

g 

I 8 = [ dy - % = = ^^(8960a 4 + 80640a 3 6 2 + 110880a 2 6 
J \Ja + by — y 2 32768 

+ 48048a6 6 + 64356 s ), 



from which we obtain 



<9 4 f , V' 2 (r)V"(r) 315tt 
dr- 



<r r 
dE 1 ' 



9E 4 J ^J2(E-V{r)) L 3 ' 
In conclusion we have 

,hs 4 f , i4r 1 1357T 1 3157T 4 71 

(-) ida A = h [——pT - „ r r , ] = h 



L 240 L 3 576 L 3 J 64L 3 
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